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ABSTRACT 



We derive an effective low energy action for QCD in 4 dimensions. The low 
energy dynamics is described by chiral fields transforming non-trivially under both 
color and flavor. We use the method of anomaly integration from the QCD action. 
The solitons of the theory have the quantum numbers of quarks. They are expected 
to be the constituent quarks of hadrons. In two dimensions our result is exact, 
namely the bosonic gauged action of WZW. 
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1. Effective action from anomaly integration 



The connection between the phenomenologically successful non-relativistic con- 
stituent-quark model (NRQM), and QCD's fundamental degrees of freedom con- 
tinues to be one of the open and puzzhng questions in particle physics. The con- 
stituent quarks which are the basic building blocks of the NRQM have the same 
quantum numbers as the much lighter and highly relativistic QCD current-quark 
fields. In addition, recent data from polarized deep inelastic lepton-nucleon scat- 
tering'^ and the ensuing conclusions about spin structure of the nucleon'*' in- 
dicate that the constituent quarks are composite objects with internal structure. 
The emergence of constituent quarks as low-energy effective excitations of QCD is 
clearly due to a nonperturbative mechanism. Finding a theoretical description of 
this mechanism is an important and interesting challenge. 

Some steps in this direction have already been made. A model combining some 
of the features of the chiral quark and skyrmion'^ approaches has been con- 
structed by Kaplan . In this model it is postulated that at distances smaller than 
the confinement scale but large enough to allow for nonperturbative phenomena 
the effective dynamics of QCD is described by chiral dynamics of a bosonic field 
which takes values in U{NcXNf). This effective theory admits classical soliton 
solutions. Assuming that they are stable and may be quantized semiclassically, 
one then finds that these solitons are extended objects with spin 1/2, and that 
they belong to the fundamental representation of color and flavor. Their mass is of 
order Aqcd and radius of order I/Aqc/j!^"' It is very tempting to identify them as 
the constituent quarks. Thus the constituent quarks in this model are "skyrmions" 
in color space. 

This is a very attractive idea, but some crucial elements of the puzzle are 
still missing, most importantly a derivation of the effective dynamics from QCD. 
However, recent work on QCD in 1+1 dimensions (QCD2) provides evidence in 
support of this picture!"' Exact nonabelian bosonization allows the Lagrangian of 
QCD2 to be re-expressed in terms of a chiral field U{x) e U{NcXNf) and of the 
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usual gauge field. This bosonized Lagrangian has topologically nontrivial static 
solutions that have the quantum numbers of baryons and mesons, constructed out 
of constituent quark and anti-quark solitons. 

When Nc solitons of this type are combined, a static, finite-energy, color singlet 
solution is formed, corresponding to a baryon. Similarly, static meson solutions are 
formed out of a soliton and an anti-soliton of different flavors. The stability of the 

mesons against annihilation is ensured by flavor conservation. These results can 
be viewed as a derivation of the constituent quark model in QCD2. Thus the idea 
of constituent quarks as solitons of a Lagrangian with colored chiral fields becomes 
exact m. D — 1-1-1. 

In this paper we will provide a constructive scheme for deriving from QCD the 
relevant effective Lagrangian. In two dimensions the procedure will yield the exact 
bosonic action of QCD2, namely the gauged WZW Lagrangian. In four dimensions 
the resulting action turns out to be a chiral Lagrangian whose leading terms are 
very similar to the action conjectured by Kaplan. 

The construction of the effective chiral action for bosonic variables is based 
on integration of the anomaly equation.'^^ The idea is to construct an action 
which will reproduce the anomaly under chiral transformations. One begins with 
the anomaly equation as the starting point and integrates it with respect to the 
chiral rotation parameter. The resulting action typically contains additional terms 
which are anomaly-free. These terms are analogous to an integration constant in 
an ordinary indefinite integral. Therefore their coefficients are arbitrary. There are 
also additional terms in the form of Chern-Simons action which are non-local. 

In order to resolve this freedom, one may adopt a strategy analogous to com- 
puting a definite integral, in this case rotating the action by a finite chiral rotation. 
This rotated action is vector invariant by choice of scheme. When the change in 
external fields is compensated by a change in the chiral phase, the rotated action 
is also axial invariant. Next, taking the difference of the rotated and unrotated 
action, one gets the effective chiral action for bosonic fields. This action includes 
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a term which is non-local, namely the Wess-Zumino term, which locally is a total 
derivative. 

We will see that in 2 dimensions this procedure generates the gauged WZW 
action, while in 4 dimensions it generates a new type of low-energy effective action 
for chiral fields transforming non-trivially under both color and flavor. We will then 
demonstrate that baryons and constituent quarks are the solitonic excitations of 
such an action. 

In order to carry out this program, we need to deal with several preliminary 
issues. First, we have to define the scheme in which the anomaly equation will 
be discussed. Having chosen the scheme, we also need to choose a suitable regu- 
larization prescription. In a given scheme the choice of the regularization cannot 
influence the results of an exact calculation. As we shall see, however, in 4 di- 
mensions the effective action contains an infinite number of terms, corresponding 
to expansion in terms of powers of external momenta divided by the cutoff. If 
we truncate the expansion at a finite number of terms and take a finite cutoff, 
the result will be regularization dependent. In 2 dimensions, however, the limit of 
infinite cutoff reproduces the exact bosonized action of WZW. 

Given an anomaly equation one can redefine the axial current and add a counter 
term which is an axial variation of a local action!^''' This creates an ambiguity in 
the definition of the quantum axial current. So we need some principles which 
remove this ambiguity. The first one is our choice of the vector scheme, we want to 
remain with a vector-conserving theory. The second one is to keep the right relation 
between the axial and the vector currents. In two dimensions the axial current is 
the dual of the vector current = —e^'^J^^. In 4 dimensions no such simple 
relation exist. We choose the heat kernel regularization!^*^' In 2 dimensions this 
procedure is well defined and the end result is a complete bosonization for gauged 
fermions. Unfortunately in 4 dimensions this cannot be done and the kinetic term 
is proportional to the cutoff squared. We therefore keep the UV cutoff finite and 
look for its physical meaning. This will be described later. The price we pay is 
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that we remain with regularization dependent terms and a low energy expansion. 
The vahie of the UV cutoff will determine the range of energies for which this 
approximation is valid. What is remarkable is that these calculations reproduce 
very accurately the numbers calculated for chiral Lagrangians in the case of flavor 

[17] 

symmetry. It is still mysterious why the heat-kernel regularization fits the results 
of experiment 

The outline of the paper is as follows. In sec. 2 we introduce the notation 
and set up the general formalism. In sec. 3 we apply the formalism to derive the 
gauged WZW action in 2 dimensions. In sec. 4 we derive the low energy effective 
action for the chiral field in 4 dimensions. In section 5 we discuss various schemes 
of choosing the group in which the chiral field is embedded, new classical solutions 
in the "[/-scheme", conserved currents including the baryon number current, and 
finally the Polyakov-Wiegmann formula extended to 4 dimensions. Its low energy 
approximation, up to dimension-4 operators, is given explicitly. 



2. General formulation 

In this section we set the notation and describe the general procedure. Those 
features which are special to 2 or 4 dimensions will be discussed in the relevant 
sections. 

The fermionic part of the QCD Lagrangian with external sources coupled to 
vector and axial currents, scalar and pseudoscalar densities labeled l^, A^, -S", P 
respectively, has the form 



Sp = - - - 5 - i^^P)^. (2.1) 

Where = V^fT", = A;^T", S = 5"T", P = P^T", Ta e G are the generators 
of some group G, normalized through Tr(r«r^) = \5°'^. In the generic case we 
think of G as the group U {Nc x Nf) , where Nc and Nf are the number of colors 
and flavors, respectively. Other groups can also be chosen, as will be discussed 
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below. V will eventually include the dynamical gauge field. For later convenience 
we introduce the chiral fields 



L 



Vu, - A 




(2.2) 



At this stage wc do not include a mass term, as this breaks the chiral symmetry 
explicitly, on top of the breaking due to anomaly. The variation of S and P, on 
the other hand, is compensated by a rotation of if) (see below). 

Since (2.1) is quadratic in the fermion fields, the fermions can be integrated 
out explicitly, yielding the determinant 



At this point we have some freedom of choice. The axial current cannot be con- 
served. We need to choose a scheme of conservation. Here, since at the end V will 

contain the dynamical field for gluons, we choose Z to be invariant under vector 
gauge transformations. 




(2.3) 



< 



5P = i [P, u] 

6L^ = d^uj + i [L^,uj] 



(2.4) 
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On the other hand, the variation of the determinant under the axial transformation 

r 5V^^i[A^,X] 



(2.5) 



5A^ = D^\ = d^\ + i[V^,X\ 
SS=- {P, A} 
5P = {S, A} 

= -9^A - i [L^, A] 
[6Ri, = di,X + t[Ri„X], 

is anomalous. This non-invariance of Z under axial transformation can be written 
in compact form as 

.SlogZ 



5A« 



[V,A,S, P] 



(2.6) 



where A is the anomaly. We integrate this equation using the chain rule 

sz 



sz _ _ sz _. 



A 



sv,, 



(2.7) 



to get an action Zi which obeys 

.SlogZi 



SX 



^ = [V,AS,P]. 



(2.8) 



Z can be written as a product Z1Z2, where Zi is not chiral invariant and can be 

calculated exactly from the anomalies, and Z2 is chiral invariant and cannot be 
calculated via this procedure. We now make a right rotation with a chiral field U, 



^ tU _ T 

rU = U-^R^U - U-^id^U 
{S + iPf = {S + iP)U 
{S - iP)^ ^ U-\S - iP) 



(2.9) 



and denote Z^ ^ Z [V^ , A^ , , P^]. Later on we will show that the vector 
invariance implies that such a rotation can be taken without loss of generality. 
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Note that U can be chosen in any subgroup of U {Nc x Nf). This point will be also 
discussed later. 

We arrive at the connection 



Sf ^ -i log Z2 



(2.10) 



^eff = -i log ^1 [V, A, S,P]+i log Zi V^,A^,S^, P' 



(2.11) 



Now let us observe the following. First, Z2 contains important physical informa- 
tion in it. It has, in principle, any power of the external sources, therefore it 
contains non-trivial information on general correlation functions. Second, Z2 does 
not depend on U. Therefore all the U physics is in SqQ. We will argue later that 
U describes the dynamics of quark solitons or baryon solitons, therefore S^q de- 
scribes their complete dynamics. Third, note that in 2 dimensions, since we get 
exact bosonization, Z2 becomes trivial. In the next sections we will treat SqQ only. 
Note that log Z^ is invariant under the transformation 



where $ represents any of the fields L^, R^, S, P in (2.9). In fact itself is 
invariant. Note however that applying the variation of $ in logZi reproduces the 
anomaly. S^q is invariant under the vector transformation 




(2.12) 



(R, 



U-^RaU - 1 



L 



U-^LM - J, 



< 



S±iP^ U-\S±iP)U 



(2.13) 



where Ju — iU dJJ. 
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We now turn to explain why the right rotation is sufficient in (2.9). The 
generahzation of (2.9) to arbitrary left and right phases is 

Rf, = 9R^Ri^9R-9R^id^9R ^^^^^ 
' {S + iPy^g2\S + iP)gR 

^{S-iPy^g^\S-iP)gL 

where g denotes {gL,9R) together. The effective action with both left and right 
rotations is 

^eff [V, A, S, P, gL, gR] = -i log [V, A, S,P]+i log Zi [R9, L^, S^, p9] . (2.15) 
Vector invariance (2.13) implies 

S,s [V, A, S, P, gL, gR] = S,s [V, A, S, P, 1, g^'gR] 
and we can define U = gJ^gR to obtain (2.11). 

3. The procedure in 2 dimensions 

In this section we calculate the effective action in 2 dimensions. In this case 
the hmit of UV cutoff going to infinity is well defined. We use Tr(757^7,y) = 2e^;y, 
eoi = 1, g'oo — ~9ii — 1- The axial variation of the action in 2 dimensions is given 

by 

F^y is given by 

F^,^. = - d^^V^ + I [V^, K] + I [A^, A^\ = ^{L^^ + R^^), (3.2) 
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and we have introduced the left and right field strengths 



LaU — dnLi, — dj^Ln + i [Ln, Lp] 

(3.3) 

R^v = dfj^Ri, — dpR^ + i [R^i Rv\ ■ 

Also, D^Ap = d^Ap + i \V^,A^. Note the term D^A^ in (3.1), which originates 
from a g^^ part rather than e^y part. This term can be calculated by recognizing 
that for ^4 = the anomaly is —■^^'^^F^^y and, for 7^ 0, F is replaced by the 
field strength of the new vector field + e^yA^ . 

Using the chain rule eq. (2.7) we have an action which has the variation (3.1) 

-I log Zi = Scs [R] - Scs [L] + ^ f MV^,Ay) - ^ [ Tr A^A^ 

2n J 2n J ^3^^ 

+ -yTr(5^) 

Where Sqs is the Chern-Simons action (an extension to 2+1 dimensions) 

Scs [^] = ^ / '"""^ T^i^^^dyA, + l^\AyA,) (3.5) 
which has the variation 

6Scs [^] = ^ / T^r{{6R^)Rp) + e^^'' Tr(i?^,(5i?,)) (3.6) 

and when using (5it!^ = d^tUJ + i [i?^, a;] the variation is 

bScs [^] = ^ / ^""^ 1V((9^i?.)u;). (3.7) 
We use the Chern-Simons anomaly contribution 



R 



u 



- Scs [R] = / Tr( J^i?,) + e^-'^ TriJ.JM {31 
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where J/j, = U ^idiJJ and = —UJiJJ ^, to calculate the effective action (2.11), 

+ ^ y {g^"" + e^'^) IV [J^i?^ - L^,U-^R^U + + L^i?^] (3.9) 
- ^ Tr I [t/-i(5 - zP)]' + [(5 + zP)[/]2 + 2(p2 _ ^2)1 

Note that it is P_ that couples to the right fermion current, and L+ to the left 
one. As can be seen the action (3.9) is the gauged WZW action of level l!^'^' We 
learn that by anomaly integration, the effective action is an exact bosonization of 
Dirac fermions. Hints for this were given in ref. [20]. Note that when P ^ 
and yS" — > M, we get a C/^ + C/~^ term, but not the expected U + term which 
corresponds to a fermionic mass term. This is because we included only effects 
that come from the anomaly, through the Jacobian of the change of variables in 
fermion sector, but not the explicit breaking coming from the mass term. 



4. The procedure in 4 dimensions 

We now turn to 4 dimensions. We follow exactly the same procedure as in the 
2 dimensional case but here we encounter some difficulties. While the result of eq. 
(3.9) is an exact bosonization of fermions coupled to vector and axial sources, its 
4-dimensional analogue is an approximation, valid only at low energies. We use 
an expansion in powers of the cutoff, assuming the higher terms are suppressed by 
powers of the cutoff squared. The scale at which the approximation breaks down 
is determined by A, the ultraviolet cutoff in the regularization procedure. 

Let us now go over to the derivation. The variation of the determinant under 
the axial transformation (2.5), including terms up to zeroth power of the cutoff, is 
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1 [21,12] 

given by 



—i- 



.5\ogZ 1 



5A« 47r2 



1 4 



+Ff,yAxAa) + -^^^^.^A^a 



+ 1 {DM^, A^A^} + 2A^D,A'^A^ + 1 [[A^, Ay] , 

+^(L>^A^A/^A'^ + A^AyD^At") + D^AyA^Af" + A^A^Df^A"^ + 0{A-^) 

(4.1) 

where H^^ = ^{R^^ — L^^). In (4.1) we have terms which are contracted with 
the metric rather than with the anti-symmetric tensor. Their coefficients are de- 
termined from a heat-kernel regularization procedure.'^^'^'^' We do not have here a 
relation like in 2 dimensions between vector and axial currents, and therefore we 
do not see simple relations among the various terms in (4.1). Using the chain rule 
(2.7), an action which has the above variation is given by 



2 a2 r 

-i log Zi = --^ y Tr (A^A^) + Sis [R] - Shs [L] 



+ 



487r2 

] 

12 



iiRnv + L^y) {Lx, Ra} + {L^Ly — -L^Ry + R^Ry)LxRa 



(4.2) 

where S^g [R] is the Chern-Simons action in 5 dimensions (see for example the 
review in Ref. [22].) 



If 3 3 

^CS [R] = 24^2 / ^'^^^'^'^ "^"^{R^idvRxdaRp + -iR^RyRxdaRp — -R^RyRxR^Rp)- 

(4.3) 

After rotation by U we get the final result eq. (2.11), 
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+ e^"^"Tr|^i?^, {Rx, J,} - R^,R,RxJa 



3 



I 

i 



y e'^'^^- Tr |i(C/-ii?^,C/ + L^,) {l^, i?^} + {L^L, - ^L^R^^ + R^R^)LxR^ 

+ j Tr [-2(D^t/-iL"^[/)2 + {D^U-^D^Uf + 2{D^D^'U-^)D^D^U 

+A{R^^Dm-^D''U + L^^Df^UD'^U-^) + U-'^Lj^^URt"' - Li^^R^""] 
+ 16 [^^,^,] + \{D^A^)^ - {A^A^f^ 

(4.4) 

where D^U = 9^C/ + iR,j,U - iULf,, Df^U'^ = df^U'^ + iLf^U''^ - iU-'^R,,. Since 
we are interested in effective action for QCD we set — R^ — in eq. (4.4), 
to obtain 



487r2 



+ 4^2 / ^''"^"Tr jiG^, [Ga, j,] - {G^G, - ^G^J, + J^J.)GaJ,} 
[ ,^.Aa i(C/-iG^,C/ + G^,) {Gx, G^} + (G^G, - ^G^G\^ + G^G^i )Gy0i 

+A{G^„,D''U-^D''U + G^.t^D^UD'^U-^) + U-^G^.r^UG'"' - G^^^G""] 

(4.5) 
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5. Bosonization Schemes and Solutions 

5.1. Different schemes 

In this section we introduce several schemes of the effective action. In the 
previous section we derived the effective low-energy action for a general coset G/H 
with H a subgroup of G. The chiral phase takes values in G, while the gauge 
field takes values in H. In QCD, H is SU{Nc). Depending on the physics one is 
interested in, there are various choices for G. We list here those which we consider 
the most important. 

1. The "[/-scheme". In this case we choose U e U{Nc x Nf). 

2. The "product scheme". The group G is SU{Nc) x SU{Nf) x U{1). One can 
parameterize U hy U = h - g, where h e SU{Nc) and g G U{Nf). 

3. The effective U{1) action. It is an effective action for the rj'. It is related to 
det 

5.2. Classical solutions in the [/-scheme 

The classical soliton solutions are embedded in an SU{2) subgroup of U {NcNf) . 
They have the form 

where r arc the generators of this SU{2) subgroup and /(r) is a radial shape 
function with boundary conditions /(O) = tt, /(oo) = 0. There are three choices 
for T. 

1. T — Ef, where Ef is a matrix in flavor space which satisfies the relation 
Ej — Ef. This is a colored soliton. This kind of soliton has been discussed 
in detail by Kaplan!^' Its baryon number is Trf(Ef). 

2. T = EcX Bf. If Ec is the identity, this is a color-singlet flavored soliton. This 
is the familiar Skyrme type soliton'*''^', widely discussed in the literature. Its 
baryon number is Trc(£^c)- 
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3. A mixed colored flavored soliton. From the general discussion in Ref. 9 and 
from the calculations in 2 dimensions'"' we expect that this new type of 
soliton corresponds to a constituent quark. 

The solitons of the theory can be taken in color space using the following 
ansatz, 

where Tc is in an SU{2) subgroup of SU{Nc) and E'y is a matrix in flavor space 
with 1 in the entry and zeroes elsewhere. In the normalization convention in 
which a nucleon has baryon number Nc, the baryon number of this soliton is 1 and 
its flavor quantum number is i, i.e. it has the quantum numbers of a quark. 

It is important to keep in mind, however, that the classical solution with the 
quantum numbers of a constituent quark has a finite energy. This is in clear 
contrast with the expectation that there are no finite energy colored states. This 
is a clear deficiency of the present treatment, and it can be traced back to the fact 
that in the anomaly integration approach, confinement has to be put in by hand, 
as it comes from the terms which we do not treat. 

One may also consider solutions which carry baryon number different from 1. 
In view of the results in 2 dimensions'"' we expect to find color-singlet solutions 
constructed out of soliton-anti-soliton pair and solutions which are constructed 
from Nc solitons. It is natural to interpret such solutions as mesons and baryons 
and their constituent solitons as constituent quarks. 

5.3. Conserved currents 

The right and left currents are given by 

if = - ^.^'"'-W, - ^AJ^r, J"} - hJ^r) (5.2) 
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and the baryon number current by the trace of their sum, 



Eqs. (5.2) and (5.3) can be conveniently written in a compact form 



SRa 



Tr 



1 



-r 



27r2 967r2 



5L,, 



= Tr 



1 



27r2 967r2 



5.4. The Polyakov-Wiegmann formula 



[24] 



Wc want to find the effective action for 2 chiral phases. The basic 
follows from the definition of the effective action SqQ (2.11). 



^eff [L, R, U1U2] = ^eff \l, R^'M + ^eff [h R, Ui] 



For zero sources we have 



S,s [0, 0, U,U2] = S,s [0, 0, Ui] + S,s [0, 0, U2] + Si [Ui, U2] 



which defines Si [Ui,U2] as 



Si [Ui, U2] = ^eff [0, - Jf , U2] - S^s [0, 0, U2] 



For 2 dimensions we have 



[24] 



•^eff [0,0, t/] = ^J Tr{J^Jn '^j e''-- Tr{J^J,J,) 
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and 



^1 = [t/l, U2] = / (^M- + ^M.) TV(J2''-^r)- (5-10) 



For 4 dimensions we have 

A2 



(5.11) 



^ i 487r^ i (5.12) 

I 3 Tjy tA to- Jl' 7A TO-\ 

+ o-^l •^2 •Jl -'2 ~ -'2 -^2 -'2 -'l J 



6. Summary and outlook 

We have derived an effective low energy action for QCD in 4 dimensions. We 
have shown that the sohtons in the theory have quantum numbers of quarks. We 
interpret these solitons as constituent quarks. An interesting plan might be to 
calculate the equations of motions numerically and compare the results to exper- 
imental values. Other quantities such as the of the constituent quarks might 
be calculated. In this case qa has a classical contribution too, in analogy with 
the Skyrme model. Multi-sohton solutions can also be considered as baryons and 
mesons constructed out of quark solitons. 
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